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ABSTRACT 

Algol (/3 Persei) is the prototypical semi-detached eclipsing binary and a hierarchical triple 
system. From 2006 to 2010 we obtained 121 high-resolution and high-S/N echelle spectra of 
this object. Spectral disentangling yields the individual spectra of all three stars, and greatly 
improved elements both the inner and outer orbits. We hnd masses of Ma = 3.39 ± 0.06 Mq, 
Mb = 0.770 ± 0.009 M© and Me = 1-58 ± 0.09 M©. The disentangled spectra also give the 
light ratios between the components in the B and V bands. Atmospheric parameters for the 
three stars are determined, including detailed elemental abundances for Algol A and Algol C. 
We hnd the following effective temperatures: Ta = 12 550 ± 120 K, Tb = 4900 ± 300 K and 
Tc = 7550 ± 250K. The projected rotational velocities are va siniA = 50.8 ± 0.8 kms“^, 
j;BsiniB = 62±2kms“^ and u© sin i© = 12.4±0.6kms“^. This is the hrst measurement of 
the rotational velocity for Algol B, and conhrms that it is synchronous with the orbital motion. 
The abundance patterns of components A and C are identical to within the measurement 
errors, and are basically solar. They can be summarised as mean metal abundances: [M/H] a = 
—0.03 ± 0.08 and [M/H]c = 0.04 ± 0.09. A carbon dehciency is conhrmed for Algol A, with 
tentative indications for a slight overabundance of nitrogen. The ratio of their abundances is 
(C/N)a = 2.0 ± 0.4, half of the solar value of (C/N)© = 4.0 ± 0.7. The new results derived 
in this study, including detailed abundances and metallicities, will enable tight constraints on 
theoretical evolutionary models for this complex system. 

Key words: stars: fundamental parameters — stars: binaries: eclipsing — stars: binaries: 
spectroscopic — stars: individual: Algol 


1 INTRODUCTION 

In his opening speach to lAU Colloquium 107 (‘Algols’), Batten 
(1989) was tempted to define an Algol system simply as “a binary 
in which the less massive component fills its Roche lobe and the 
other, which does not, is not degenerate”. The resulting evolution¬ 
ary paradox can be solved by postulating an episode of mass trans¬ 
fer from what was initially the more massive component to its then 
less massive, and consequently less evolved, companion (Crawford 
1955). This first, rapid, phase of mass exchange between the com¬ 
ponents eventually leads to an Algol-type system. Numerical cal¬ 
culations of the evolution of stars in binary systems have proved 
this hypothesis to be very plausible (cf. Paezynski 1971). 

Direct comparison of the observed properties of Algol sys¬ 
tems to evolutionary model calculations are sparse, because only a 
few Algols have reliable measurements of their physical properties 
(Maxted & Hilditch 1996), and because models are unavailable for 
the relevant parameters. The main observational obstacle in the de¬ 
termination of the accurate stellar quantities for Algols arises due 
to the intrinsic faintness of the Roche-lobe filling mass-losing com¬ 


ponent, the reason they were detected only due to technological 
advances in light detectors (c.f. Tomkin 1989). Model calculations 
for Algol with both conservative and non-conservative binary evo¬ 
lution have only been calculated by Sama (1993). 

Changes in the chemical composition of the surface layers of 
stars which have experienced mass transfer are a feature of theoret¬ 
ical predictions. Layers that were once deep inside the star become 
exposed after mass transfer. This makes it possible to observation- 
ally probe CNO nucleosynthetic processes in stellar cores, and the 
efficiency of different mixing processes in stellar interiors. In that 
sense CNO abundances can serve as a sensitive probe of stellar 
structure and evolution. 

Carbon depletions have been detected in Algol and in Algol- 
type systems (Parthasarathy, Lambert & Tomkin 1979, 1983; 
Cugier & Hardorp 1988; Cugier 1989; Tomkin et al. 1993, Ibanoglu 
et al. 2012). However, comparisons to theoretical models were 
rather scarce since detailed evolutionary calculations were per¬ 
formed for only a few Algol-type systems, and were limited to 
tracing the evolution of carbon (De Greve & Cugier 1989; Sama 
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1992, 1993; De Greve 1993). On the observational side, Algol sys¬ 
tems can be rather difficult. The secondary stars are Roche-lobe 
filling subgiants which are usually intrinsically faint, but still bright 
enough to contaminate the spectra of the primary stars. Many con¬ 
tain third components, as is the case for Algol itself, which fur¬ 
ther complicates spectroscopic analyses. Therefore, we initiated a 
new observational project to obtain high-resolution echelle spec¬ 
troscopy of Algol systems to which we could apply the powerful 
technique of spectral disentangling (hereafter SPD) and thus obtain 
the individual spectra of the component stars. This makes detailed 
spectroscopic analyses much more straightforward, and allows the 
study of the spectral characteristics of all components over wide 
wavelength ranges. To aid the comparison with theoretical evo¬ 
lutionary models of semi-detached binary systems, we have also 
considered other diagnostic species. The primary target is nitrogen 
abundance, which is very sensitive to CNO nucleosynthetic pro¬ 
cesses. Theoretical models of binary star evolution have also been 
calculated for non-conservative mass transfer. A detailed descrip¬ 
tion of our methodology was given in a study of the hot Algol u Her 
(Kolbas et al. 2014). a large sample of mass transferring binaries 
of Algol type can reveal all the different evolutionary paths which 
these complex binary systems encounter. 

The chemical composition was determined from spectra cov¬ 
ering the full optical range for only a few Algol-type binaries. Be¬ 
side uHer (Kolbas et al. 2014), abundance patterns have been de¬ 
rived for RZ Cas (Narusawa et al. 2006, Tkachenko et al. 2009), 
TW Dra (Tkachenko et al. 2010), TX UMa (Glazunova et al. 2011), 
and AS Eri (Narusawa 2013). All these objects except u Her have a 
pulsating primary component. 

In this work we present new high-quality optical spectroscopy 
of Algol with two main goals Firstly, improved measurement of the 
masses of all three components through spectral disentangling, with 
a special emphasis in studying the spectrum of a faint Roche-lobe- 
filling component over a broad spectral range. Secondly, spectro¬ 
scopic determination of the effective temperature, abundance pat¬ 
terns and bulk metallicities from disentangled spectra of the indi¬ 
vidual components. Our aim is to assemble useful observational 
constraints for comparison to evolutionary model calculations; in 
particular, the first measurement of the bulk metallicity of the com¬ 
ponents of Algol. 

We first summarise the rich observational history of Algol, and 
highlight some recent observational results. In Sect. 3 we present 
new echelle spectra of Algol. Sect. 4 covers the SPD of the spectra 
and the determination of new orbital elements for both inner and 
outer orbits. The analysis of the disentangled spectra for all three 
components is presented in Sect. 5, including detailed calculations 
of abundances for Algol A and metallicity for Algol C. The results 
are discussed in the context of mass transfer between Algol A and 
B in Sect. 6. Finally, Sect. 7 summarises our results and discusses 
future work. 


2 ALGOL IN A NUTSHELL 

Algol (/3Persei, HD 19356) is a hierarchical triple star system 
(Frieboes-Conde et al. 1970; Sdderhjelm 1980). It is extremely 
bright {V = 2.12 mag) and its periodic variability is the subject 
of a long observational history. There is evidence that Algol’s peri¬ 
odicity was recorded by the ancient Egyptians three millennia ago, 
in the Cairo Calendar (Jetsu et al. 2013). The first well-documented 
discovery of the 2.867 day periodicity, and an explanation in terms 
of stellar eclipses, was given by Goodricke (1783) based on naked- 


eye observations. Thanks to its brightness, Algol also has a long 
history of observations with many techniques and at many wave¬ 
lengths. An outstanding and exhaustive account of these studies 
was given by Wecht (2006). 

The inner pair consists of a late-B-type star in orbit with an 
early-K-type subgiant which fills its Roche lobe. This close sys¬ 
tem exhibits partial eclipses, and is the prototype of the Algol class 
of eclipsing binaries. A tertiary component pursues a 680 d or¬ 
bit around the inner pair; its spectral classification has been vari¬ 
ously given as late-A, early-F, and Am. This description of Algol 
emerged from decades of primarily spectroscopic and photometric 
studies (cf. Friebos-Conde et al. 1970; Hill et al. 1971; Sdderhjelm 
1980; Richards et al. 1988). The main difficulty of these early stud¬ 
ies was the intrinsic faintness of the cool subgiant, which is exac¬ 
erbated by dilution due to the light of component C. Algol C is 
brighter than Algol B, but its contribution has previously been very 
difficult to quantify (Richards et al. 1988). A breakthrough came 
with its eventual spectroscopic detection in the infrared (Glusheva 
& Esipov 1967), and then radial velocity (RV) measurements which 
led to the determination of the dynamical masses for the compo¬ 
nents (Tomkin & Lambert 1978). 

Other observational techniques have been helpful in constrain¬ 
ing the orbital and physical characteristics of Algol. Rudy & Kemp 
(1978) found phase-locked polarization and independently deter¬ 
mined the inclination of the eclipsing binary. Subsequently Kemp 
et al. (1983) discovered eclipse polarization, or the Chandrasekhar 
effect, further evaluated by Wilson & Liou (1993). The problem 
of the determination of the mutual orientations of the orbits in the 
Algol system has remainined open (Kemp et al. 1981). 

Since the successful detection of Algol C with speckle inter¬ 
ferometry (Labeyrie et al. 1974), Algol has often been a target for 
interferometric measurements, which have recently settled the is¬ 
sue of the true orientation of both orbits. Csizmadia et al. (2009) 
spatially resolved the inner pair with long baseline interferometry 
in combination with very long baseline interferometry (VLBI) ra¬ 
dio observations. They found the inner orbit to be prograde, in dis¬ 
agreement with the retrograde movement found by Lestrade et al. 
(1993) from radio observations. This long-term controversy was 
solved by Zavala et al. (2010), who simultaneously resolved all 
three stellar components in the optical, achieving the then highest 
precision in angular measurements. The outer orbit was shown to 
be prograde and the inner orbit retrograde. VLBA radio (Peterson 
et al. 2010, 2011) and CHARA iT-band measurements (Baron et 
al. 2012) subsequently confirmed these orientations for both orbits. 
The most recent highlight of the interferometric studies is an un¬ 
ambiguous spatial resolution of the three stars in the Algol system 
with angular resolution <0.5 mas (Baron et al. 2012). This made 
possible the determination of the orbital and physical character¬ 
istics of the components (angular sizes and mass ratios) indepen¬ 
dently from previous studies. Baron et al. (2012) also determined 
the mutual inclination of the orbits to be much closer to perpen¬ 
dicularity than previously established. The distortion of the Roche- 
lobe-filling component B is clearly seen in the reconstructed image. 

Another highlight of recent work on Algol is the first three- 
dimensional reconstruction of this system (Richards et al. 2012). 
These three-dimensional tomograms revealed previously unde¬ 
tected evidence of the mass transfer process, such as loop promi¬ 
nences and coronal mass ejections. Early predictions of the super¬ 
hump phenomenon in Algol, i.e. the gas between the stars in close 
pair being threaded with a magnetic field even though the hot mass- 
gaining star is not know to have a magnetic field (Retter et al. 2005) 
have been supported by this new technique. Algol B is a late-type 
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Figure 1. Depiction of the phase distribution of the observed spectra for 
the inner (upper panel) and outer (lower panel) systems. The RV curves 
(solid lines) were calculated from the orbital elements found via SPD (Ta¬ 
ble 5). Filled circles represent the RV expected for each observed spectrum, 
and are not measured RVs. This is because the SPD analysis returns the or¬ 
bital elements directly, without passing through the intermediate step of RV 
determination (cf. Simon & Sturm 1994). Spectra obtained during eclipse 
are not shown due to the rotational distortion seen in the primary eclipse 
(cf. Sect. 4). Phases for the inner eclipsing system were calculated with the 
ephemeris of Kreiner (2004) with the primary minimum at phase 0. For the 
outer system the ephemeris from Baron et al. (2012) was used, where Tq is 
the time of periastron passage. 

magnetically active subgiant found to be a strong radio and X-ray 
source (Wade & Hjellming 1972; White et al. 1986; Stern et al. 
1992, 1995), and dominates in the spectrum of Algol at these wave¬ 
lengths. In particular, X-ray studies have made possible a determi¬ 
nation of the chemical abundances in the corona of Algol B (Drake 
2003) as well both the X-ray-bright star Algol B and the X-ray-faint 
star Algol A (Yang et al. 2011). 


3 HIGH-RESOLUTION SPECTROSCOPY 

Our observational programme of high-resolution and high signal- 
to-noise (S/N) spectroscopy of Algol was initiated with two observ¬ 
ing runs in 2006 and 2007 at the 2.5-m Nordic Optical Telescope 
(NOT) at La Palma, Spain. We obtained 85 spectra of Algol us¬ 
ing the Fibre-fed Echelle Spectrograph (FIES; Telting et al. 2014). 
FIES is housed in a dedicated climate-controlled building and has 
a high thermal and mechanical stability. The wavelength scale was 
established from thorium-argon exposures taken regularly through¬ 
out the observing nights. We used fibre 4 in bundle B, giving com¬ 
plete spectral coverage in the interval 3640-7360 A at a reciprocal 
dispersion ranging from 0.023 A px“^ in the blue to 0.045 A px“^ 


in the red. The resolution of the instrument is roughly 3.5 px, giving 
a resolving power of 48 000. An exposure time of 15 s was used for 
all spectra, resulting in continuum S/N ratios in the region of 200- 
500 in the B and V bands. 

In order to cover the long-period outer orbit of Algol, spec¬ 
troscopic observations were taken from 2009 to 2010 at the Bo- 
hyunsan Optical Astronomy Observatory (BOAO), South Korea. 
A set of 36 spectra were secured with Bohyunsan Optical Echelle 
Spectrograph (BOES) mounted on the 1.8-m telescope (Kim et al. 
2007). BOES has multiple spectral resolving powers up to 80000, 
and covers thee wavelength range 3600 to 10 200 A. A thorium- 
argon lamp was used for wavelength calibration and S/N ratios of 
300-550 were achieved. 

The spectra were bias-subtracted, flat-fielded and extracted 
with the IRAiQ echelle package routines. Normalisation and merg¬ 
ing of the orders was performed with great care, using custom pro¬ 
grams, to ensure that these steps did not cause any systematic errors 
in the resulting spectra (see Appendix). 


4 ORBITS AND MASSES THROUGH SPECTRAL 
DISENTANGLING 

The method of SPD enables determination of the orbital elements of 
binary and multiple systems, along with the simultaneous separa¬ 
tion of the individual spectra of the components, in a self-consistent 
way (Simon & Sturm 1994). In that sense it is a generalisation 
of the Doppler tomography method of Bagnuolo & Gies (1991) 
which was the first successful reconstruction of individual spectra 
for binary stars but which relies on predetermined RVs for the com¬ 
ponents for each spectrum. In SPD only a time series of observed 
spectra are needed, with roughly uniform coverage of the RV mo¬ 
tion of the stars (cf. Hensberge et al. 2008). Optimisation of the RVs 
is bypassed in favour of directly fitting for the orbital elements, as 
first implemented by Simon & Sturm (1994). 

With no need for input template spectra, SPD does not suffer 
from biases due to template mismatch (cf. Hensberge & Pavlovski 
2007 and references therein). However, no comprehensive study of 
the error propagation in SPD and comparison to cross-correlation 
has yet been undertaken. Some initial studies have indicated the 
method is ‘well-behaved’ (Hynes & Maxted 1998; Ilijic et al. 2001; 
Hensberge & Pavlovski 2007; Southworth & Clausen 2007). 

The SPD of Algol is challenging for the following reasons: (i) 
the system is triple; (ii) the secondary component is barely discem- 
able in the optical spectra; (iii) the period of the outer orbit is long 
(680 d); and (iv) eclipse spectra could not be used due to significant 
distortion of the line profiles due to the Rossiter-McLaughlin effect 
(Rossiter 1924; McLaughlin 1924). 

Spectra of several multiple systems have previously been suc¬ 
cessfully separated using different variants of the SPD method (cf. 
Fremat et al. 2005; Gonzales et al. 2006; Lee et al. 2008; Tamajo et 
al. 2012). Faint components with a fractional light contribution of 
only ~5% have previously been detected using SPD (e.g. Pavlovski 
et al. 2009; Lehmann et al. 2013; Tkachenko et al. 2014; Borkovits 
et al. 2014), including cases with only ~2.5% (Torres et al. 2014), 
and 1.5-2% (Holmgren et al. 1999; Mayer et al. 2013) in the V 
spectral region. As was shown in our preliminary reports on SPD 
of Algol (Pavlovski et al. 2010; Kolbas et al. 2012) we are pushing 

^ IRAF is distributed by the National Optical Astronomy Observatory, 
which ai'e operated by the Association of the Universities for Research in 
Astronomy, Inc., under cooperative agreement with the NSF. 
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Figure 2. Two portions of disentangled spectra for the three components of the Algol system. The plots are not on the same scale, and disentangled spectra for 
Algol B and C are shifted up by 0.04 and 0.10 of the continuum flux, respectively, for clarity. 


the limit of the method because of the extreme faintness of Algol B 
(see Sect. 5.1). 

The rotational distortion of the spectral lines during eclipse 
(the Rossiter-McLaughlin effect) violates the principal assumption 
of spectral disentangling that spectral line profiles are not intrinsi¬ 
cally variable, so spectra obtained during eclipse cannot be used. 
Without such spectra there is insufficient variation in the light con¬ 
tributions of the components during the orbit, so the zeroth mode 
in the Fourier expansion is undetermined. This has the consequence 
that spectra can be only separated, and remain in the common con¬ 


tinuum of the system. Renormalisation of the separated spectra of 
the individual components is discussed further in Sect. 5.1. 

To perform SPD we used the code FDBINARvQ (Ilijic et al. 
2004). FDBINARY implements disentangling in the Fourier do¬ 
main (Fladrava 1995) using fast fourier transforms, making it 
computationally efficient. It can also account for two orbits and 
three individual spectra, as is required for Algol. Five spectral 
segments were selected, paying particular attention to including 


^ hhtp://sail.zpf.fer.hr/fdbinary 
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Table 1. The final solution for the orbital elements of the inner and outer 
orbit in the Algol triple system as obtained by SPD. The periods are taken 
from Baron et al. (2012), as well as the eccentricity of the outer orbit. Trial 
calculations for the inner orbit has shown it is circular, and in the final cal¬ 
culations eA—B and uja—b were set to 0. 


Quantity 

Notation 

Unit 

Value 

Error 

Inner orbit 

Period 

Pa-b 

d 

2.871362 

fixed 

Eccenti'icity 

eA-B 


0. 

fixed 

Periastron longitude 

^A-B 

deg 

90. 

fixed 

RV semiamplitude 

Ka 

km s“^ 

44.1 

0.2 

RV semiamplitude 

Kb 

km s~^ 

194.2 

1.2 

Mass ratio 

Qa-b 


0.227 

0.005 

Outer orbit 

Period 

Pab-c 

d 

680.168 

fixed 

Time of periastron 

7ab-c,o 

d 

2454433.2 

1.1 

Eccenti'icity 

eAB-C 


0.227 

fixed 

Periastron longitude 

WAB-C 

deg 

138.1 

0.6 

RV semiamplitude 

Kab 

km s~^ 

11.9 

0.4 

RV semiamplitude 

Kc 

km s~^ 

32.9 

0.8 

Mass ratio 

Qab-c 


0.364 

0.009 


enough spectral lines of Algol B to allow its velocity semiampli¬ 
tude to be precisely determined. For the blue spectral region 
this is challenging since the fractional light contribution is less 
than 1% (see Sect. 5.1 and Fig. 2). The following spectral seg¬ 
ments were used: 4176.9-4277.0, 4400.2^614.7, 4929.1-5118.9, 
5140.6-5305.1 and 5311.0-5520.0 A. Only spectra taken outside 
eclipse were used, giving a total of 112 input spectra. To limit the 
computational demand we also stacked spectra taken in close suc¬ 
cession, resulting in a total of 49 input spectra. Such spectra typi¬ 
cally covered a time interval of 3 min, for which the RV of Algol 
B changes by a maximum of 0.5kms“^, so velocity smearing is 
not a problem. The FIES and BOES datasets were also analysed 
first separately with the outer orbit fixed according to recent inter¬ 
ferometric results (Czismadia et al. 2009; Zavala et al. 2010; Baron 
et al. 2012). The solution for the inner orbit was then used as an 
initial set of the orbital elements for solving both orbits simultane¬ 
ously using the 49 spectra. 

Hill et al. (1971) found that the inner orbit is slightly eccentric, 
ba-b = 0.015 ± 0.008. In our SPD calculations the eccentricity of 
the inner orbit always converged to ca-b = 0.0 with a high confi¬ 
dence. Therefore, in the final set of the calculations, the eccentricity 
of the inner orbit was fixed to zero. 

Whilst the coverage of the orbital cycle for the inner orbit 
is good (Fig. 1), the phase distribution of the observations for the 
outer orbit is not sufficient to determine all parameters simultane¬ 
ously, so we fixed the eccentricity of this orbit. Recent independent 
interferometric studies give a very consistent result on its value of 
cab-c = 0.227±0.002 (Baron etal. 2012). Since the simplex al¬ 
gorithm is used for minimisation in FDBINARY, care must be taken 
to avoid becoming trapped in local minima. We used 50 runs with 
1000 iterations each. Fixing bab-c made the convergence stable 
and consistent. Our final solution for the two orbits is given in Ta¬ 
ble 1. It represents the mean values from the solutions of SPD in the 
five selected spectral segments specified above. The errors quoted 
in Table 1 are standard deviations of the mean calculated from the 
solutions obtained for the five segments. More sophisticated error 
calculations for SPD, such as the jackknife method (Pavlovski et al. 
2011, 2014), were beyond our computational resources. 


The comparison between the results presented in Table 1 and 
those from previous studies is interesting. First, we did not confirm 
the eccentricity for the inner orbit (ba-b) which was reported in 
some previous studies, but rejected in others. In an extensive spec¬ 
troscopic study Hill et al. (1971) examined a series of medium- 
resolution photographic spectra obtained at the Dominion Astro- 
physical Observatory. They found ba-b = 0.015±0.008 but were 
not able to detect the spectrum of the secondary star. It is not clear if 
and how this might affect the determination of RVs for the primary. 
Also, it should be noted that Hill et al. (1971) attributed a 32 year 
periodicity to apsidal motion in the inner orbit. The possibility of 
orbital eccentricity for the eclipsing pair has not been supported by 
past photometry of Algol (Wilson et al. 1972; Guinan et al. 1976; 
Sdderhjelm 1980; Kim 1989). 

The velocity semi-amplitude for Algol A found by Hill et al. 
(1971), K\ = 44.0 ± 0.4kms“^, agrees within the uncertainties 
with the value derived in our work, 44.1 ± 0.2kms“^. Also, it is 
encouraging that their velocity semi-amplitudes of the outer orbit, 
Kab ~ 12.0 ± O.dkms”*^ and Kc = 31.6 ± 1.2kms“^, are in 
agreement with our values, TTab = 11-9 ± 0.4kms“^ and Kc = 
32.9 ± 0.8kms“^. Hill et al. (1971) were also able to determine 
the eccentricity of the outer orbit, bab-c = 0.23 ± 0.04, which 
was corroborated and improved by interferometric measurements 
(Zavala et al. 2010; Baron et al. 2012). 

Two portions of disentangled spectra of the three components 
are presented in Fig. 2. The secondary spectrum is clearly isolated 
in the whole optical range, even though it contributes only 1% of 
the light around 4500 A and somewhat less than 2% around 5500 A 
(Sect. 5.1). The RV semi-amplitude for Algol B, Kb, has two al¬ 
most equally deep local minima around 194 and 211 kms“^. It is 
difficult to trace the reasons for this ambiguity but, besides the 
small light contribution, these might include imperfect phase dis¬ 
tribution of the observations for the inner orbit, gaps in the phase 
distribution for the outer orbit, and possible systematics due to the 
use of two spectrographs and varying observing conditions. The 
for the lower solution is slightly better, so we adopted this as our 
final solution, Kb = 194.2 ± 1.2kms“^. This value agrees with 
the only previous measurement, Kb = 201 ± 6km s“^ (Tomkin & 
Lambert 1978). 

Tomkin & Lambert (1978) calculated the spectroscopic mass 
ratio using the Ka determined by Hill et al. (1971) and their own 
Kb measurement, finding qa-b = 0.219 ± 0.007. This mass ratio 
was subsequently used in many studies on Algol, as it was the first 
mass ratio found via dynamical effects. Our value is slightly higher, 
qA-B = 0.227 ± 0.002. The mass ratio can be estimated from 
light curve for semi-detached eclipsing binaries, and there have 
been attempts to derive it this way. Kim (1989) used the ‘g-search 
method’ to obtain exactly to the same value from BV light curves, 
gphot = 0.227±0.002. Determination of the astrometric orbit from 
the interferometric measurements also can yield the mass ratio. 
Baron et al. (2012) solved the orbital elements for both orbits and 
the mass ratio of the inner orbit from the relative positions of the 
components derived from the visibility functions. In addition, they 
derived the mass ratio of the outer orbit independently from paral¬ 
lax assumptions. They found qa-b = Mb/Ma = 0.219 ± 0.017 
and qAB—c ~ Me/{M a ± Mb) = 0.456 ± 0.022. This qa-b 
agrees with our value to within its relatively large uncertainties. 

The mass ratio for the outer orbit derived from these different 
techniques is highly discrepant and needs clarification. However, it 
is encouraging that spectroscopic analyses are consistent, as Hill et 
al. (1971) found qab-c = 0.380 ± 0.051, compared to our value 
of qAB—c = 0.362 ± 0.042. 


© 2014 RAS, MNRAS 000,[T]-?? 






6 V. Kolbas et al. 



Figure 3. Determination of the fractional light contribution for Algol A from optimal fits of He I lines in its disentangled spectrum. Synthetic NLTE He I line 
profiles were calculated for Tefj = 12 550 K, log g = 4.05 and v sini = 50.5 km s“^. 


Table 2. Summary of the fractional light contribution of the components of 
Algol to the total light of the system determined by different techniques. 
‘Photometry’ indicates light curve analysis. Only B and V photometry are 
listed as this is the wavelength region covered by our spectta. References: 
(1) Wilson et al. (1972); (2) Demircan (1978); (3) Richards et al. (1988); 
(4) Kim (1989); (5) Zavala et al. (2010); (6) This work. 


Method 

B 

V 

Source 

Photometry 

0.940 0.010 0.050 

0.893 0.029 0.078 

(1) 

Photometry 

0.909 0.041 0.050 

0.857 0.065 0.078 

(2) 

Photometry 

0.929 0.012 0.059 

0.894 0.035 0.071 

(3) 

Photometry 

0.854 0.040 0.106 

0.814 0.067 0.119 

(4) 

Interferometry 


0.863 0.063 0.065 

(5) 

Spectroscopy 

0.943 0.008 0.049 

0.915 0.018 0.067 

(6) 


5 ATMOSP HE RIC DIAGNOSTICS 

5.1 Renormalisation of disentangled spectra 

Without eclipse spectra or significant variations in the fractional 
light contribution between the components, SPD can only be per¬ 
formed in the separation mode. The renormalisation of the indi¬ 
vidual spectra of the components can then be done using external 
information on the individual components’ fractional contribution 
to the total light of the system. 

A simple way is to use the light ratio between the components 
from the light curve solution for eclipsing binaries (cf. Hensberge 
et al. 2000). However, the light curve solution of Algol is rather 
uncertain because the eclipses are partial and there is ‘third light’ 
from component C. On top of these obstacles, it is extremely diffi¬ 
cult to acquire photometry of Algol due to its brightness and lack 
of suitable comparison stars within 10°. 

It is therefore not surprising that only a few complete light 
curves for Algol have been published. The most recent photometry 
(Kim 1989) was a series of photoelectric measurements in the B 
and V filters, obtained from 1982 to 1984 at Yonsei Observatory, 
Korea. The other light curves of Algol from the era of photoelec¬ 
tric photometry are those of Guinan et al. (1989) and Wilson et al. 
(1972). In addition, Stebbins secured an incomplete light curve in 
multiple passbands from 1949 to 1951 (Stebbins & Gordon 1975), 
Al-Naimiy et al. (1985) obtained unpublished observations in 1981, 
and light curves in the ultraviolet (Eaton 1975), and near-infrared 
(Chen & Reuning 1966) exist. 

Previous light curve analyses have not yet given a unique so¬ 


lution for Algol, primarily because of its complexity. Solutions of 
partially eclipsing binary stars inevitably suffer from degeneracy, 
which should be lifted using additional information (Garcia et al. 
2014; Southworth et al. 2007). This is illustrated in Table 2 in which 
fractional light contributions (light dilution factors) for Algol are 
listed from several works. Wilson et al. (1972) used an estimate 
of the third light contribution derived spectroscopically by Fletcher 
(1964). The same approach was followed by Demircan (1978) who 
analysed the photometry published in Wilson et al. (1972). Since 
the secondary component had not yet been detected in the spec¬ 
trum of Algol, the mass ratio between the components of the inner 
eclipsing pair had to be assumed from other considerations. But it 
is evident from the solutions derived in Richards et al. (1988) and 
Kim (1989) that the dynamical mass ratio measured by Tomkin & 
Lambert (1978) did not help to provide a consistent description of 
the Algol system. Richards et al. (1988) calculated a grid of models 
with different assumptions for the radiative properties of the third 
star, whilst Kim (1989) attempted to derive its contribution directly 
from his BV light curves. In this context we also list estimates 
of the light ratios measured in the interferometric observations by 
Zavala et al. (2010). Their estimates suffered from large uncertain¬ 
ties, and the fractional light contribution of ~6% for Algol B is 
certainly too large. 

Facing all these uncertainties in the determination of the frac¬ 
tional light contributions of the components, we rely on the extrac¬ 
tion of these quantities from the disentangled spectra themselves. 
The information on the light dilution factor is contained in the spec¬ 
tral line depths in the disentangled (separated) spectra of an indi¬ 
vidual component. But caution is needed at this point since other 
effects could change the line depth in stellar spectrum and mimic 
light dilution effects in binary or multiple systems, such as metal 
abundance or chemical peculiarity of the component(s). Algol C is 
often characterised as an Am star, beginning with the spectroscopic 
analysis by Fletcher (1964). 

In the optical spectrum of Algol the most prominent lines are 
HI and He i originating in the photosphere of the primary compo¬ 
nent (cf. Struve & Sahade 1957). Since Hei lines are present in 
only the primary star’s spectrum, we used them for the determina¬ 
tion of its effective temperature. Teg, a, and fractional light con¬ 
tribution, If A- We used the program STARFIT (Kolbas et al. 2014) 
which compares disentangled component spectra to a grid of calcu¬ 
lated theoretical spectra. STARFIT can handle the following param¬ 
eters: effective temperature Tea, surface gravity logg, fractional 
light contribution If, projected rotational velocity tisini, relative 
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velocity shift between disentangled spectrum and rest-frame labo¬ 
ratory wavelengths, vo, and continuum corrections, cc. Optimisa¬ 
tion is performed by a genetic algorithm (Charbonneau 1995), and 
can be done in constrained mode with disentangled spectra of the 
components for which the sum of the fractional light contributions 
should be 1 (Tamajo et al. 2011). Grids of theoretical spectra 
were calculated in LTE (local thermodynamic equilibrium) using 
the UCLSYN code (Smith 1992, Smalley et al. 2001) 

Tefj^A is constrained to be in the range 12000 to 13 000 K 
(Richards et al. 1988; Kim 1989). Our disentangled spectra con¬ 
tain three Hei lines suitable for the optimisation, at 4471, 4713 
and 5015 A. Unfortunately, the Hei lines at 4388, 4920, 5047, and 
5788 A are either on the wings of Balmer lines, or are contami¬ 
nated by metal or telluric lines, so are not suitable. To facilitate 
optimisation we first determined v sin i by optimal fitting of un¬ 
blended metal lines in the primary’s disentangled spectrum (v sin i 
or FWHM does not depend on the light dilution). In hot stars, a de¬ 
generacy exists between Teff and logg for Hi and Hei lines, so 
we fixed logg in our calculations to loggA = 4.05 (Richards et 
al. 1988). Also, it should be noted that model atmospheres are cal¬ 
culated assuming the ‘standard’ (solar) helium abundance with the 
fraction of helium atoms A^He/(AlH + A^He) = 0.089. 

We found that the three He i line profiles are not well and con¬ 
sistently reproduced in the LTE approximation. This is not a sur¬ 
prise because discrepancies in the Hei line profiles are the most 
severe at the low and high end of the temperature sequence (Auer 
& Mihalas 1970, 1973). Thus we decided to calculate Hei line 
profiles in non-LTE. Our grid of theoretical spectra (Kolbas et 
al. 2014) was extended down to Tefj = 12 000 K. This is be¬ 
low the usual borderline used for NLTE spectrum synthesis of 
Teff — 15 000 K. Almost perfect fits were achieved with Hei 
line profiles in NLTE (Fig. 3, solid lines). Optimisation was per¬ 
formed separately for each line to find the wavelength-dependent 
fractional light contribution. An excellent convergence was found 
for almost same effective temperature. Teg ,a = 12 600 ± 90 K, 
(/(4471-f4713) = 0.943±0.002, and 1/(5015) = 0.915±0.002. 
Uncertanties in Teff and If were calculated with a Markov chain 
Monte Carlo technique (Ivezic et al. 2014). A refined Teff mea¬ 
surement for Algol A is discussed below (Sect. 5.2). 

As described in Sect. 4, SPD enabled us to reconstruct almost 
the entire optical spectrum of Algol B, the first time this has been 
achieved (Figs. 2 and 4). A first attempt to match it to theoreti¬ 
cal spectra with reff,B ~ 4500 K indicated a low fractional light 
contribution, notably less than indicated by light curve solutions 
(Table 2). Again, the program STARFIT was used with only sur¬ 
face gravity as a fixed parameter, log(/B = 3.11 (Richards et al. 
1988). The projected rotational velocity converged to vb sintB = 
62 ± 2kms“^, in good agreement with the expected synchronous 
rotational velocity for this star (Waynch = 61.8 ± 0.5kms”*^. The 
optimal fractional light contributions of Algol B are: Idf (4500) = 
0.008 ± 0.001 and ld/(5500) = 0.018 ± 0.001, with Teg.B = 
4900 ± 300 K. This is the most unexpected result of our study but 
quite secure since the spectral lines for the late G or early K sub¬ 
giant are intrinsically deep and, for the expected v sin i, the light 
dilution is certainly no more than 2% in V and less than 1% in B. 

5.2 Effective temperatures for the components 

Detailed spectral analysis is possible for Algol A and Algol C, mak¬ 
ing it possible to fine-tune their Teff measurements. The values de¬ 
termined in Sect. 5.1 served as starting points. A spectral analysis 
of the entire disentangled spectral range (4400-5800 A) was per- 





Wavelength [A] 

Figure 4. Optimal fits of the disentangled spectra of the components in Al¬ 
gol. From top to bottom are shown Algol A, Algol B, and Algol C, respec¬ 
tively. These short portions of the components’ disentangled spectta show 
the quality of the optimal fitting for the atmospheric parameters listed in 
Table 5, except for Algol B (middle panel) for which best fits for two TgffS 
are shown (solid line 5000 K, and dotted line 4000 K). 


formed, using UCLSYN. The equivalent widths (EWs) were mea¬ 
sured for suitable lines in the renormalised disentangled spectra. 
Model atmospheres were calculated in LTE with the program AT- 
LAS9 (Kurucz 1979), and abundances calculated for the measured 
EWs. Since Fe lines are the most numerous in the spectra of both 
stars, they served for the determination of Teff and microturbulent 
velocity, ^t. The spectrum of Algol A contains many Fell lines, 
but those of Fe i are few and weak. That of Algol C contains many 
Fe I lines but also a substantial number of Fe ii lines, which means 
fine-tuning of Teff ,c is possible using the Fe ionisation balance (see 
Fig. 4). 

We determined Teff and in a few iteration steps. Teff ,a was 
tuned so there was no correlation between Fe ii abundance and the 
excitation potential, EP, whilst the microturbulence, ^t,A was found 
by requiring the Fell abundance to be independent of EW. We 
found Teff,A = 12 550 ± 120 K and Ct,A = 0.4 ± 0.2 km s"^ The 
uncertainties were calculated from the uncertainties in the determi¬ 
nation of the slopes of log e(Fe ii) versus EP for TeS and log e(Fe ii) 
versus EW for ^t. 

The atmospheric parameters for Algol C were determined 
in the same manner. Moreover, the Fe ionisation balance could 
be used for fine-tuning Teff.c- Our final results are Teff^c = 
7 540 ± 80K, and ^t,c = L64 ± 0.08kms“^. The uncertain¬ 
ties were derived in the same way as for Algol A, except that 
Fe I was used. The ionisation balance of Fe is well satisfied with 
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Figure 5. Dependence of iron abundance on wavelength for Algol A (filled 
red circles) and Algol C (filled blue circles). In the case of Algol A, the 
abundances of Fe II are shown. For Algol C, Fe I abundances are shown. The 
scatter is large for Algol C because the S/N of the renormalised disentangled 
spectrum of this component is much less than that for Algol A due to the 
differing fractional contributions to the total light of the system. 
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Figure 6. Abundance pattern measured for the photospheric composition 
of Algol A (symbols) compared to the standard solar composition (grey 
shading) which represents the Icr uncertainty in the solai' composition de¬ 
termined by Asplund et al. (2009). 


a difference in abundance derived from Fe i and Fe ii lines of only 
Aloge(Fe) = 0.05 ±0.11. 

Since we have renormalised the disentangled spectra of the 
components with the light dilution factors determined from these 
spectra themselves, the wavelength dependence of the Fe abun¬ 
dance could be an important check of the correctness of the pro¬ 
cedure. Fe II and Fe i lines are well distributed in the optical spectra 
of Algol A and Algol C, respectively. This test is more sensitive for 
Algol C as the multiplication factor needed to normalise its spec¬ 
trum to unit continuum is much higher than in the case of Algol 
A, -^Ib.S versus -^1.1. As is illustrated in Fig. 5 no wavelength 
dependence of iron abundance is present for either of the compo¬ 
nents. This provides encouraging support for the reliability of our 
estimates for the luminosity contribution of the components to the 
total light of the system. The for Algol A is also in agreement 
with that generally found for late B-type stars (Fossati et al. 2009), 
which is usually considerably less than 1 km s~^. 


5.3 The elemental composition and metallicity 

The number of spectral lines available for the determination of 
the photospheric chemical composition of Algol A are rather lim¬ 
ited. This does have the advantage that line blending and smearing 
do not prevent detailed analysis despite the relatively high v sin i 
of 50kms“^. Beside the strong hydrogen (not disentangled) and 
helium lines, the most prominent lines are of Fe ii. For the Tefr 
of Algol A, spectral lines of C, N and O are rather weak; the 
EWs of C lines are all <4.5 mA except for the Cii 4267 A line. 


Table 3. Photospheric abundances derived for Algol A. Abundances aie 
expressed relative to the abundance of hydrogen (column 2), loge(i7) = 
12.0. The third column gives the number of lines used, and abundances 
relative to the Sun are given in the fourth column. 


El. 

log e(X) 

N 

[X/H] 

C 

8.64±0.02 

3 

0.21±0.05 

N 

7.92±0.04 

4 

0.09±0.06 

0 

8.93±0.02 

11 

0.24T0.05 

Mg 

7.49±0.03 

7 

-0.10±0.05 

Si 

7.46±0.04 

14 

-0.05±0.05 

P 

5.38±0.05 

3 

-0.03±0.05 

S 

7.02±0.02 

21 

-0.10±0.06 

Ca 

6.37±0.02 

3 

0.05±0.04 

Ti 

4.91±0.03 

11 

-0.02±0.09 

Cr 

5.62±0.04 

9 

0.00±0.07 

Mn 

5.48±0.05 

5 

0.06±0.06 

Fe 

7.39±0.02 

70 

-0.08±0.05 

Ni 

6.18T0.02 

17 

-0.02±0.05 


Table 4. Photospheric abundances derived for Algol C. Explanations are 
the same as for Table[3] 


El. loge(X) 


c 

8.53±0.04 

0 

8.72±0.06 

Na 

6.31 ±0.08 

Mg 

7.50±0.09 

Si 

7.58±0.06 

S 

7.16±0.04 

Ca 

6.37±0.09 

Sc 

3.03±0.07 

Ti 

5.01±0.07 

V 

4.10±0.10 

Cr 

5.69±0.07 

Mn 

5.33±0.08 

Fe 

7.45±0.08 

Co 

5.03±0.07 

Ni 

6.25±0.07 

Cu 

4.19±0.13 

Y 

2.32±0.06 

Zr 

2.59±0.06 

Ba 

2.38±0.34 

La 

1.35±0.07 

Ce 

1.69±0.06 

Nd 

1.57±0.08 


N [X/H] 


21 

0.10±0.06 

6 

0.03±0.08 

7 

0.10±0.09 

7 

-0.09±0.10 

23 

0.07±0.07 

10 

0.04±0.05 

21 

0.02±0.09 

12 

-0.13±0.08 

106 

0.08±0.08 

17 

0.21±0.13 

97 

0.07±0.08 

23 

-0.09±0.09 

296 

-0.02±0.09 

21 

0.10±0.09 

161 

0.05±0.08 

3 

0.01±0.14 

14 

0.11±0.08 

11 

0.09±0.07 

3 

0.29±0.18 

6 

0.23±0.08 

16 

0.11 ±0.07 

8 

0.15±0.09 


<3 mA for N, and <15 mA for O. Therefore, abundance estimates 
for these elements should be considered with caution, except for 
C II 4267 A whose EW of 44.8 mA gives a carbon abundance of 
loge(C) = 8.27 ± 0.06. This is 0.16 ± 0.08 dex less than so¬ 
lar (Asplund et al. 2009). A deficiency of carbon is also indicated 
from the comparison of the EW for the C ii 4267 A line to EWs of 
standards from Tomkin et al. (1993) and Ibanoglu et al. (2012). 

The mean photospheric metal abundance of Algol A has been 
calculated for all elements listed in Table 3 except CNO. The abun¬ 
dances of these elements were not used because of the possibility 
of changes caused by mass transfer, and the small EWs of the lines 
on which they are based. The remaining ten elements give an av¬ 
erage abundance of [M/FIJa = —0.03 ± 0.08. This is less than 
la away from the [Fe/H]A = —0.08 ± 0.05 value derived from 
the most numerous lines. We calculated the bulk metallicity using 
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Figure 7. Abundance pattern measured for the photospheric composition of 
Algol C (black dots) compared to the average abundance pattern of ‘normal’ 
A-type stars (upper panel, open triangles) and Am stars (lo\ver panel, open 
squares), following Gebran et al. (2010). 


the approximate equation, Z = Zq 10 , finding it to be solar 

to within its uncertainty: Za = 0.014 ± 0.002 (using the present 
day solar metallicity Zq = 0.0134 from Asplund et al. 2009). The 
measured elemental abundances in the photosphere of Algol A are 
displayed in Fig. 6 relative to standard solar values (Asplund et al. 
2009). 

The rich and relatively unbroadened (vc sin ic ~ 12 km s“^) 
spectrum of Algol C made it possible to determine the photospheric 
abundances for 20 elements. Seven appear in two ionisation stages 
(neutral and singly-ionised). In the spectral range studied (4500- 
5700 A) the most numerous lines are of Fe i, Ni i, Cr i and Ti ii. The 
measured abundances are listed in Table 4, along with the num¬ 
ber of spectral lines used, and abundances relative to standard solar 
(Asplund et al. 2009). The slightly larger uncertainties compared to 
Algol A, despite having more lines available, are due to the smaller 
S/N of the renormalised disentangled spectrum after multiplication 
by ~16 to put the continuum level to unity. In Fig. 7 we compare 
the photospheric composition of Algol C with the abundance pat¬ 
tern derived in Gebran et al. (2010) for a sample of ‘normal’ A-type 
stars (upper panel) and Am stars (lower panel). The abundances in 
Algol C do not show the principal characteristics of the Am phe¬ 
nomenon, i.e. strong Sc underabundance, often a Ca underabun¬ 
dance, and a moderate to strong overabundance of iron-group ele¬ 
ments. The mean metallicity for Algol C is consistent with solar, 
at [M/H]c = 0.04 ± 0.09. The corresponding bulk metallicity is 
Zc = 0.012 ± 0.002. 

The photospheric elemental compositions of Algol A and Al¬ 
gol C are the same to within the la uncertainties, with a formal 
difference of A[M/H] = [M/H]^ - [M/HJg = -0.07 ± 0.12. 
Both correspond to the standard solar composition (Asplund et al. 
2009). 


6 DISCUSSION 

6.1 Fundamental properties of the components of Algol 

With the velocity semiamplitudes determined for all the compo¬ 
nents of the Algol system, we are able to derive their dynami¬ 
cal masses. The orbital inclinations are taken from Richards et al. 
(1988) for the inner orbit, (a-b = 81.4 ± 0.2, and from Baron 
et al. (2012) for the outer orbit, (ab-c = 83.66 ± 0.03. In the 
latter study ia-b was also derived, and was in agreement with 
Richards et al. (1988) but with lower accuracy. These inclinations 
are well constrained from several different studies and techniques. 
We also used the orbital periods from Baron et al. (2012): Pa-b = 
2.867328 ± 0.00005 d and Pab-c = 680.168 ± 0.54 d. We find 
the masses Ma = 3.39±0.06 Mq, Mb = 0.770±0.009 M© and 
Me = 1.58 ± 0.09 Mq. From the outer orbit we get the sum of 
the masses for Algol A and Algol B, Mab = (Ma ± MB)outer = 
4.38 ± 0.27 Mq, which is in accordance with the sum of the in¬ 
dividual masses of Algol A and Algol B from the inner orbit, 
(Ma ± MB)inner = 4.16 ± 0.06 Mq. Table 5 contains the phys¬ 
ical properties of all three stars in the Algol system derived in the 
current work. 

The component masses we find are about 8% smaller than the 
most commonly quoted values (Richards et al. 1988). This corrob¬ 
orates the findings of Baron et al. (2012), although these authors 
found masses about 15% small than those from Richards et al. 
(1988). Besides improving the mass values we were able to im¬ 
prove the measurement precision to 1.8% for Algol A and 1.2% for 
Algol B. That for Algol C is less improved, at 5%. 

Published determinations of the radii of the stars suffer from 
degeneracy due to the third light and partial eclipses (see Sect. 
5.1). Richards et al. (1988) found Ra = 2.90 ± 0.04 Rq, Rb = 
3.5 ± 0.1 Rq and Rc ~ 1.7 Rq (no error given). Baron et 
al. (2012) achieved better than 0.5 mas spatial resolution in the 
7/-band with the CHARA interferometer, and unambiguously re¬ 
solved the three stars. They found angular diameters of (f)A = 
0.88±0.05 mas, = 1.12±0.07 mas and <(>c = 0.56±0.10 mas. 
With the parallax of Algol determined by Zavala et al. (2010), tt = 
34.7 ± 0.6 mas, they found linear radii of Ra = 2.73 ± 0.20 Rq, 
Rb ~ 3.48 ± 0.28 Rq and Rc = 1.73 ± 0.33 Rq. The interfer¬ 
ometric measurements agree with our results to within their large 
uncertainties. 

The dimensions of the Roche-lobe filling component B are 
constrained by the mass ratio. This gives another way of determin¬ 
ing its linear radius and synchronous rotational velocity. With our 
mass ratio for the semidetached pair, qa-b = 0.227 ± 0.005 we 
have a relative radius of Algol B of tb = 0.251 ± 0.002 (ra¬ 
dius at the sides of the secondary star perpendicular to the line 
connecting the centres of the two stars). The semi-major axis of 
the inner orbit is oa-b = 13.65 ± 0.07 Rq. The linear radius 
is Rb = tb fflA-B = 3.43 ± 0.01 Rq. This is consistent with 
the Baron et al. (2012) value. The synchronous rotational veloc¬ 
ity for this Rb is Uaynch.B = 60.2 ± 0.2kms“^. We measured 
ub sin iB = 62±2 km s“^ from the spectral line broadening, which 
also supports our estimate of the radius of Algol B. 

If we use the condition of synchronous rotation for Algol 
A, which does not have to be fulfilled, then for the measured 
UAsintA = 50.5 ± 0.8kms“'^ we get Ra = 2.87 ± 0.04Rq. 
If Algol A has been spun up by mass transfer (cf. Packet 1981; 
Dechamps et al. 2013) then this would be an upper limit for its ra¬ 
dius. An improvement in direct interferometric measurements of its 
angular diameter, and/or revising the light curve analysis with new 
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Table 5. Astrophysical quantities for the components of Algol triple system derived in this work (masses, Tgfjs, nsinj values, or calculated in conjuction 
with previous studies (radii, logg values, synchronous velocities). The radius for Algol C is from the interferometric study hy Baron et al. (2012). All other 
quantites are determined or calculated in this work. 


Quantity 

Notation 

Unit 

Algol A 

Algol B 

Algol C 

Mass 

M 

Mq 

3.39 ±0.06 

0.770 ± 0.009 

1.58 ±0.09 

Radius 

R 

Rq 

2.87 ±0.04 

3.43 ±0.01 

1.7 ±0.3 

Surface gravity 

logs 

[cgs] 

4.05 ±0.01 

3.254 ± 0.006 

4.18 ±0.16 

Effective temperature 

Teff 

K 

12 550 ±120 

4900 ± 300 

7550 ± 250 

Projected rotational velocity 

V sin i 

km s~^ 

50.5 ±0.8 

62 ± 2 

12.4 ±0.6 

Synchronous velocity 

^synch 

km s~^ 

51.2 ±0.5 

64 ± 1 

14.1 ± 2.5 

Microturbulence velocity 


km s“^ 

0.4 ±0.2 


1.68 ±0.06 

Mean metal abundance 

[M/H] 


-0.03 ±0.08 


0.04 ± 0.09 

Bulk metallicity 

Z 


0.014 ±0.002 


0.012 ± 0.002 


constraints from the spectroscopy presented in this work, might al¬ 
low the radius of this component to he measured to high precision. 

6.2 Chemical composition and evolution of the components 

Predictions of carbon depletion in the atmospheres of mass¬ 
transferring systems prompted several observational studies. Car¬ 
bon underabundances relative to solar were found by Cugier & 
Hardorp (1988) in an analysis of far-UV spectra secured from 
the lUE satellite. The carbon abundance they reported for Al¬ 
gol A, relative to the modern standard solar value, is [C/H] = 
—0.32 ± 0.20 dex. Tomkin et al. (1993) studied a sample of Algol- 
type binaries, including Algol itself, and confirmed carbon de¬ 
ficiencies in the whole sample. For Algol A they found [C/H] 
= —0.22 ±0.15 dex using high-resolution observations of the 
CII 4267 A line. The reported value is relative to the average abun¬ 
dance, loge(C) = 8.28 ± 0.21, they determined for the standard 
stars. Our result reported in Sect. 5.3, [C/H]a = —0.16 ± 0.08, 
confirms a small carbon depletion for Algol A, both in terms of the 
abundance determined from the strongest carbon line in its spec¬ 
trum (CII 4267 A) and the deviation of its EW from a calibration 
derived from standard late-B type stars. Nitrogen lines show a slight 
overabundance with log €{N) = 7.97 ± 0.02, but with the caution 
that the N ii lines used are very weak (EW = 1.0-2.7 mA). Refer¬ 
ring to the standard solar nitrogen abundance (Asplund et al. 2009), 
the abundance deviation for Algol A is [N/H]a = 0.14 ± 0.05. 

The ratio between C and N abundances is a sensitive indica¬ 
tor of CNO nucleosynthesis and the efficiency of mass-transfer and 
mixing processes in stellar interiors. For Algol A we get (C/N)a = 
2.0 ± 0.3 which, compared to the solar value, (C/N)© = 4.0 ± 0.7, 
indicates that a change in the C/N ratio in Algol A has been de¬ 
tected. CNO-processed layers from Algol B, formerly the more 
massive of the two inner components, are evidently now exposed on 
the surface of Algol A, the mass-gaining and currently more mas¬ 
sive star. In our previous study on the hot Algol-type binary sys¬ 
tem uHer (Kolbas et al. 2014) we determined C/N = 0.89 for the 
mass-gaining component, u Her contains components with masses 
of 7.8 and 2.8 Mq, substantially larger than the masses of the inner 
pair in Algol. As predicted by detailed chemical evolution mod¬ 
els, Kolbas et al. (2014) found a stronger carbon depletion in the 
mass-losing component than in the mass-gaining companion. The 
progenitor mass of what is now the less massive component in 
uHer was estimated from evolutionary model calculations to be 
MdonoT 7.2 Mq, whilst an estimate for the initial mass of Algol 
B is MB.init = 2.7 Mq (Sarna 1993). 

Whilst we succeeded in separating and reconstructing the 


spectrum of Algol B, detailed analysis is premature. The disentan¬ 
gled spectrum still suffers from a low S/N, and the spectral lines 
have high rotational broadening. However, this object dominates 
the X-ray spectrum of the system. Drake (2003) used Chandra Low 
Energy Transmission Spectrograph observations to determine the 
abundances of C and N in the corona of Algol B. The analysis was 
performed relative to the ‘standard’ star HR 1099. These two stars 
have shown many similarities in their X-ray spectra, except for the 
strengths of C and N lines. For N an enhancement by a factor of 
three compared to the standard star is found, whilst no C lines are 
detected in the Algol B spectrum, indicating a C depletion relative 
to ‘standard’ by a factor of ten or more. Moreover, Drake (2003) 
found a standard (solar) Fe abundance for Algol B, in fine agree¬ 
ment with our findings for Algol A and Algol B. 


7 CONCLUSION 

The first systematic echelle spectroscopic survey of Algol, a hier¬ 
archical triple system, has been conducted from 2006 to 2010 and 
covers the entire optical spectral range. The technique of spectral 
disentangling makes it possible, for the first time, to reveal the indi¬ 
vidual spectra of all three components in Algol. In our SPD analysis 
the orbital elements for both orbits were optimised, leading to im¬ 
proved measurements of the masses of all three component stars. 
The uncertainties for the masses of the inner (eclipsing) system are 
now below 2%. 

A detailed spectroscopic analysis of the disentangled spectra 
of the individual components was undertaken, and yielded an accu¬ 
rate determination of the effective temperatures and projected rota¬ 
tional velocities for all three stars, and the photospheric elemental 
abundances and bulk metallicity for Algol A and Algol C. Equally 
importantly, tight constraints on the fractional light cotribution of 
all three components have been derived. The spectroscopic light ra¬ 
tio for the inner partially-eclipsing binary could be used to lift the 
degeneracy in the determination of the radii from light curves of 
the system. This would make the system an excellent test case for 
theoretical work on modelling the evolution of mass-transferring 
binaries, resulting in an improved comprehension of Algol systems. 
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APPENDIX A: BLAZE FUNCTION CORRECTION 

Extracted and wavelength-calibrated echelle spectra (Sect. 2) are 
split into orders, each of which has a response function determined 
by the blaze shape. In each order, the intensity of signal smoothly 
rises from the edges to the centre of frame. It is due to the optical 
setup and characteristics of echelle spectrograph, and has sine form 
(see Barker 1984). Blaze function removal from extracted spectral 
orders is a critical step in the reduction of echelle spectra. Any sys¬ 
tematic errors in the blaze function correction produces periodic 
ripples in the shape of the continuum, which in turn affects the 
depth and thus equivalent widths of spectral lines. A major obsta¬ 
cle to removal of the blaze function is the presence of broad spectral 
lines, in particular Balmer lines in late-B and early-A stars. These 
can be wider than a single spectral order, so determining the blaze 
function using points on the spectral continuum is impossible. 

The most widespread method to correct for the blaze func¬ 
tion is to divide the spectrum with a blaze function obtained from 
flat-fields (c.f. Erspamer & North 2002; Skoda et al. 2008, and ref¬ 
erences therein). However, this often gives a wavy shape in the 
merged spectral orders, which then still have to be locally nor¬ 
malised. For unknown reasons blaze functions of stellar and flat- 
field spectra are not identical, and undesired ripples appear and de¬ 
grade the quality of merged spectra. 

One of the most important prerequisites for spectral disentan¬ 
gling is the quality of the normalised merged echelle spectra. To 
avoid systematic errors, we remove blaze functions in more thor¬ 
ough and interactive way. First, the spectral orders suitable for fit¬ 
ting of the blaze function are selected (in Fig. Al these are orders 
#39, and 42). These are orders without broad spectral lines and with 
a substantial amount of continuum. High-degree polynomial func¬ 
tions are fitted by least-squares through selected continuum points. 
Wavelengths of maximum intensity of the blaze functions are also 
calculated, and used for the interpolation of the positions of blaze 
function head in the remaining orders. Blaze functions are now in¬ 
terpolated to these orders, as illustrated the upper panel of Fig. Al. 
The results of the correction using interpolated blaze functions are 
shown in the middle panel of Fig. Al. The bottom panel of Fig. Al 
illustrates the quality of interpolation and merging between two se¬ 
quential orders containing a Balmer line. It should be emphasized 
that overlapping parts fall on the edges of spectral orders, and suf¬ 
fer a lower S/N than central parts. Despite this, merging is good to 
much better than 1%. 
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Wavelength [a] 


Figure Al. Blaze function correction and merging the spectral orders con- 
taning broad Balmer lines. From the top to the bottom: Identification of 
suitable orders for the determination of the blaze function, and interpolation 
to the positions of echelle orders containing a Balmer line (upper panel), 
merging of spectral orders after removing htted (orders #39 and #42) and 
interpolated (orders #40, and #41) blaze functions (middle panel), and close 
up of merged orders in the wings and core of a Balmer line (bottom panel). 
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